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1. Introduction 

Classical and quantum information processing network architectures utilize light (optical 
photons) for the transmission of information over extended distances, ranging from 
hundreds of meters to hundreds of kilometers [TJ |2] . The utility of optical photons stems 
from their weak interaction with the environment, large bandwidth of transmission, 
and resiliency to thermal noise due to their high frequency (~ 200 THz). Acoustic 
excitations (phonons), though limited in terms of bandwidth and their ability to transmit 
information farther than a few millimeters, can be delayed and stored for significantly 
longer times and can interact resonantly with RF-microwave electronic systems [3J. 
This complimentary nature of photons and phonons suggests hybrid phononic-photonic 
systems as a fruitful avenue of research, where a new class of optomechanical circuitry 
could be made to perform a range of tasks out of reach of purely photonic and phononic 
systems. A building block of such a hybrid architecture would be elements coherently 
interfacing optical and acoustic circuits. The optomechanical translator we propose in 
this paper acts as a chip-scale transparent, coherent interface between phonons and 
photons and fulfills a key requirement in such a program. 

In the quantum realm, systems involving optical, superconducting, spin or charge 
qubits coupled to mechanical degrees of freedom [U El [TJ El |9] have been explored. 
The recent demonstration of coherent coupling between a superconducting qubit and a 
mechanical resonance by O'Connell, et al. [10], has provided an experimental backing 
for this vision and is the latest testament to the versatility of mechanics as a connecting 
element in hybrid quantum systems. In the specific case of phonon-photon state transfer, 
systems involving trapped atoms, ions, nanospheres [TT1 [T2T [T3 1 [T H 1TB] . and mechanically 
compliant optical cavity structures [TS] have all been considered. In these past studies, 
the state of an incoming light field is usually mapped onto the motional state of an atom, 
ion, or macroscopic mirror, through an exact timing of control pulses, turning on and off 
the interaction between the light and mechanical motion in a precise way. The ability to 
simultaneously implement phononic and photonic waveguides in optomechanical crystal 
(OMC) structures [TTj opens up the opportunity to implement a traveling-wave phonon- 
photon translator (PPT). Such a device, operating continuously, connects acoustic and 
optical waves to each other in a symmetric manner, and allows for on-the-fly conversion 
between phonons and photons without having to precisely time the information and 
control pulses. In effect, the problem of engineering control pulses is converted into a 
problem of engineering coupling rates. 

Our proposal for a PPT is motivated strongly by recent work [TSJ [Tj5] on radiation 
pressure effects in micro- and nano-scale mechanical systems [201 ED, E21 [231 EU EBJ EE] . 
Furthermore, the concrete realization of a PPT is aided by the considerable advances 
made in the last decade in the theory, design and engineering of thin-film artificial quasi- 
2D (patterned membrane) crystal structures containing photonic [271 EH E3 EH EU [32] 
and phononic (33J [Ml ESJ [351 EZ] "band gaps" . Such systems promise unprecedented 
control over photons and phonons, and have been separately subject to extensive 
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investigation. Their unification, in the form of OMCs which possesses a simultaneous 
phononic and photonic bandgap [TTJ |38l [39j SQj SI] , and in which the interaction between 
the photons and phonons can be controlled, promises to further expand the capabilities 
of both photonic and phononic architectures and forms the basis of the proposed PPT 
implementation. 

The outline of this paper is as follows. In Sections [2] and [3] we introduce and study 
the PPT system as an abstraction, at first classically and then quantum mechanically. 
After introducing the basic system, its properties and its scattering matrix, we study 
the effects of quantum and classical noise on device operation. In Section [| we design 
and simulate a possible physical implementation of the system, utilizing recent results 
in simultaneous phononic-photonic bandgap materials [39]. Finally, in Section [5j we 
demonstrate a few possible applications of the PPT. Focusing first on "classical" 
applications, we evaluate the performance of the PPT when used for the implementation 
of an optical delay line and wavelength converter. Finally, we show how such a system 
could be used in theory to do high fidelity quantum state transfer between optical and 
superconducting qubits. 

2. Outline of Proposed System 




Figure 1. Full system diagram. Circles represent resonant modes, while rectangles 
represent waveguides. Blue is for photonics, beige is for phononics, a color-scheme 
followed in other parts of the paper. The coupling h between the two optical modes is 
modulated by the intervening phonon resonance. 

The proposed PPT system, shown in Fig. [TJ consists of a localized mechanical 
resonance (b) which couples the two optical resonances (01,02) of an optomechanical 
cavity via radiation pressure. External coupling to and from the mechanical resonance 
is provided by an acoustic waveguide, while each of the optical resonances are coupled to 
via separate optical waveguides. Multi-optical-mode optomechanical systems have been 
proposed and experimentally studied previously in the context of enhancing quantum 
back-action, reduced lasing threshold, and parametric instabilities (HI SSI SU SSI SSI 
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H7J HE]- Here we use a two-moded optical cavity as it allows for the spatial filtering and 
separation of signal and pump optical beams while reducing the quantum noise in the 
system, as is explained below. 

A general description of the radiation-pressure-coupling of the mechanical and 
optical degrees of freedom in such a structure is as follows. For each of the two high-Q 
optical resonances of the cavity we associate an annihilation operator dk and a frequency 
Uk (k = 1,2). Geometric deformation of the optomechanical cavity parameterized by 
x, changes the frequencies of the optical modes by gk{x). The deformation, due to the 
localized mechanical resonance with annihilation operator b and frequency Q, can be 
quantized and given by x = xzpp(& + W). There is also a coupling between the two 
optical cavity modes given by h(x), where for resonant intermodal mechanical coupling 
the cavity structure must be engineered such that Q = uj 2 — u±. In a traveling wave PPT 
device consisting of the two optical cavity resonances and a single mechanical resonance, 
the lower frequency cavity mode (ai in this case) is used as a 'pump' cavity which enables 
the inter-conversion of phonons in the mechanical resonance (b) to photons in the second, 
higher frequency, optical cavity mode (02) through a two-photon process in which pump 
photons are either absorbed or emitted as needed. The 'signals' representing the phonon 
and photon quanta to be exchanged will thus be contained in b and 01, respectively. 

As described, the Hamiltonian of this system is, 

H = huj\d\ai + hw2a\a2 + MWb + h(gi(x)a\di + g2(x)ala 2 ) 

+ nh(x)(alax + a\a 2 ) + i- s /2K he E pump (e~ iu ' Lt a\ + e lij3Lt a x ), (1) 

where we have added a classical optical pumping term of electric field amplitude E pump 
and frequency ul- Optical pumping is performed through one of the optical waveguides 
with (field) coupling rate to the a% cavity resonance given by Ki >e . In addition to the 
waveguide loading of each optical resonance (nk,e), the total optical loss rate of each 
cavity mode includes an intrinsic component of field decay due to radiation, 

scattering, and absorption. Similarly for the mechanical resonance, we have a field 
decay rate given by 7 = 7 e + 7« which is a combination of waveguide loading and intrinsic 
losses. The constant parts of both h(x) and gk{x) (MO), 9k(fy) can be eliminated by a 
change of basis and are thus taken to be zero. As discussed below, it is advantageous to 
choose a cavity structure symmetry in which g^(x) = up to linear order in x. In fact, 
we can generally assume that the mechanical displacements are small enough to make 
the linear order the only important term in the interaction. Assuming then a properly 
chosen cavity symmetry, 

g k (x) = g ■ (b + P) = and h(x) = h ■ (b + 6 j ), 

which yields a simplified Hamiltonian 

H = fajjia\a\ + hbJ2a\a2 + Klb% + hh(b + tf)(alai + d\a2) 

+ z V / 2^^ P ump(e-^*a t 1 + e^'ai). (2) 

Treating the system classically and approximately, we can write each intracavity 
photon and phonon amplitude, and their inputs (see Fig. [2]) as a Fourier decomposition 
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Figure 2. Optical sidebands and processes for phonon-photon translation. The optical 
pump is located on sideband ai.o, on resonance with the first cavity mode at frequency 
There are three relevant optical sidebands to consider in the translation process, 
ot.2 _|_, ai.o and a 2 ,-- The inter-sideband photon scattering gives rise to phonon emission 
and absorption. The state transfer occurs through scattering between sidebands a 2:+ 
and a^O) whereas inter-sideband scattering between a\fl and a 2 - can be thought of 
as phonon noise. Note that for the 5 = case considered here there are no sidebands 
at uj\ + ft for cavity mode a\ and no sidebands at uj\ for mode a 2 . 



of a few relevant sidebands: 

ai(t) = a lfl e- iullt + a^+e^ 1 ^ 1 + a^e^ 1 ^ (3) 

a 2 {t) = a 2fi e- iullt + c^+e^" 1 ^* + aa^e^ 1 ^ (4) 

b(t) =p + P + e- iSlt (5) 

b- m (t) = ft n , + e- iQt (6) 

a 2 , in (t) = « in , + e-^ 1+Q ) t (7) 

b mt {t) = f3 out , + e- tnt (8) 

a 2 ,out(t) = « out , e-^* + a out , + e-^ 1+ ^ + a^e^^. (9) 



The equations of motion arrived at from the system Hamiltonian (presented generally 
in the following section) then become algebraic relations between the a and (3 sideband 
amplitudes. By ignoring the self-coupling term (g = 0), pumping on-resonance with 
cavity mode a\ (ljl = Ui), and engineering the optical cavity mode splitting for 
mechanical resonance (u 2 = u± + Q) , we arrive at classical sideband amplitudes, 

a lj+ = a X - = a 2 , = (3 = 0, (10) 
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From here we see that the central sideband amplitude of cavity mode a\, a^o, is 
proportional to the sum of a term containing the pump field -E pum p and terms containing 
products of the optical and mechanical sideband amplitudes. By increasing E pump the 
effect of the other sidebands on the pump resonance amplitude can be made negligible, 
and we assume here and elsewhere in this work that the pump sideband is generally left 
unaffected by the dynamics of the rest of the system. As desired the optical sideband 
which contains mechanical information is a2,+ since it is the only sideband directly 
proportional to /3+. The constant of proportionality between these two terms is seen to 
contain both h and a^o, demonstrating the role of the pump beam in the conversion 
process. Since coherent information transfer between the optics and mechanics is 
occurring between {3 + and a 2 ,+ , it is desirable to remove the effects of the lower energy 
photonic sideband, ot2,-- This sideband can be made significantly smaller in magnitude 
than ai2,+ in the sideband-resolved regime where Q » k. A convenient way to visualize 
all of the processes in the system is shown in Figure [2] where the photonic sideband 
amplitudes «2,+ and a\ t Q are represented as "energy levels", with transitions between 
them being due to the emission and absorption of phonons. 

From this approximate analysis it is clearly suggested that in a sideband resolved 
optomechanical system, a state-transfer process is possible between the phononic and 
photonic resonances, and the process is controlled by a pump beam [16]. A more in-depth 
study of the system dynamics required to understand how such processes may be used 
for traveling wave phonon-photon conversion, and a full investigation of all relevant 
noise sources required to understand the applicability of such a system to quantum 
information, is carried out in the following section. 

3. Analysis 

A detailed treatment of the operation of a traveling phonon-photon translator is carried 
out in this section. At first, the dynamics of the system are simplified while still 
taking into account the noise processes related to the sideband «2,-. In this way 
one is left with an effective 'beam-splitter' Hamiltonian, which describes the coherent 
interaction between the optics and mechanics, while the aforementioned noise processes 
are accounted for through an effective increase in the thermal bath temperature. This 
is followed by a treatment of the traveling-wave problem through a scattering matrix 
formulation, which provides insight into the role of the intracavity pump photon number 
(|ai,o| 2 ) and optimizing the state transfer efficiency. 
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3.1. Simplified Dynamics of the System 

Starting from the Hamiltonian in equation ([2]), a set of Heisenberg-Langevin Equations 
can be written down, [33]: 



a i 
6 

The input coupling terms as written above include both external waveguide coupling 
and intrinsic coupling due to lossy channels (see Fig. [I]). Separated, the intrinsic (with 
subscript i) and extrinsic (with no subscript) components look as follows, 



V27^ n = V^fkn + V^^in.i, (19) 

K k = K>k,i + K k,e, (20) 

7 = li + 7e- (21) 

As the fluctuations in the fields are of primary interest, each Heisenberg operator can 
be rewritten as a fluctuation term around a steady-state value, 

ai(t) — (ai + a 1 )e^ lt , (22) 

a 2 (t) -> (« 2 + a 2 )e-^ 2 *, (23) 

6(t) — (/3 + 6)e~ in *. (24) 

Assuming that the pump beam is driven resonantly with a\ (ujl = wi), the c-number 
steady-state values are equal to (ax,at2,/3) = _E pump , 0, 0). 



For the fluctuation dynamics, with A = (w 2 — Wi) — tt, the resulting equations are: 



pump 1 



(ki + zwi)ai — z/i(6 + 6 T )a 2 + ^JIkx^E, 
(k 2 + 1002)02 — ih(b + b^)ai — V2K 2 a' 2 in , 

i06 — ih(a) 2 a i + &1&2) — 76 + 76^ — ^/2/yb[ 



^ai )in ,(15) 
(16) 

(17) 



fix = - «ioi - iha 2 (be- t2nt + S^e^* - V^a^e^ 1 ' (25) 
a 2 = - K2 a 2 - i/i( ai + ai)(S + 6t e + i2nt ) e + iA * _ ^a 2iin e iW2 * (26) 

S = - 7 (S - Pe +imt ) - ih(at + a 1 )a\e + ^ +2 ^ t 

- ih{a x + aO^ae^ - ^b[ n e int (27) 

Ignoring all the mechanically anti-resonant terms for now, and invoking the rotating 
wave approximation (RWA) valid when A « Q and Q » \acx\h, we arrive at the simplified 
set of fluctuation equations, 

ai = - Kiax - iha 2 Pe~ iAt - V2^a[ - n e iu>lt (28) 
h 2 = - k 2 cl 2 - ih(&! + ai)6e +iA * - V2na' 2 - m e iuJ2t (29) 

b = - 76 - ih(&i + aOW^' - ^fb' in e im (30) 
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By ignoring all of the counter-rotating terms proportional to e +l2Qt , we have also 
neglected the noise processes alluded to previously due to the a 2 - sideband. Of the 
mechanically anti-resonant terms which have been dropped, the terms proportional to ct\ 
(hai& 'e +l2nt in equation (26) and haia\e +%mt in equation (27)) are the largest and most 
significant in terms of error in the rotating wave approximation. These terms correspond 
to aiaffi + h.c. in the Hamiltonian and cause inter-sideband photon scattering between 
the pump, ai o, and its lower frequency sideband, a 2 - as shown in Figure [2] This inter- 
sideband scattering process causes emission and absorption of phonons in the mechanical 
part of the PPT, thus in principle, even when the extrinsic phonon inputs are in the 
vacuum state, spontaneous scattering of photons from ai >0 to a 2 ,- m ay populate the 
mechanical cavity with a phonon. 

This effect was studied in Refs. [HOI EI] in the context of quantum limits to 
optomechanical cooling. Similar to that work, a master equation for the phononic 
mode with the u>i — Q optical sideband adiabatically eliminated results in an additional 
phononic spontaneous emission term given by, 

Pb,s P on = —7^ (2Ppb - Wp - pW) (31) 

K (f ) + 1 

where G = h\af\. The master equation for the mechanics, found by tracing over all 
bath and optical variables, is of the form p = pt,,i + Pb,s P on + Pb,c — i/h~[Hb, p], where the 
terms on the right hand side of the equation are respectively the intrinsic phononic loss, 
the phonon spontaneous emission, the phonon-waveguide coupling, and the coherent 
evolution of the system. The first two terms can be lumped together into an effective 
intrinsic loss, p' bi = p b ,i + p b , sp0 n, 

p' b . = 7 '.(n' + \){2bpP - b r bp - pb r b) + 7i 'n'(26 t p& ~ bP p - pbP), (32) 

li = Tt(! ~ ™spon), (33) 
-/ Tl -\- Ti spon /oa\ 

n = — ' ( 34 ) 

-L '''spon 

and 

_ G 2 1 

nSP ° n = ^(f) 2 + l' 
Hence, by assuming that the intrinsic loss phonon bath is at a modified temperature with 
occupation number n', and changing the intrinsic phonon loss rate to 7 4 ', the spontaneous 
emission and intrinsic loss noise are lumped into one effective thermal noise Liouvillian 
for the mechanics. Note that it is possible in this model to have 7^ negative when 
^s P on > 1; however, the motional decoherence rate is always positive and given by 
ifi = li{n + n spon ). 

In Section |3.2| we will see that the optimal state transfer efficiency is given by 
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G 2 = 7« [f] in which case 

nspon * j { feJ ' (36) 
in the sideband resolved limit. In the case where g # (recall this is the self-coupling 
radiation pressure term), photons scattering from the pump into the lower frequency 
sideband (u\ — O) can scatter into the a\ cavity mode which is only detuned by Q. This 
is to be compared with the g = case considered thus far in which the detuning is 2Q 
for scattering into the ai mode. As such, for the g # case there will be roughly four 
times the spontaneous emission noise, with nfj^° % 'y/-yi(K/fl) 2 . 

A final simplification can be made by neglecting the fluctuations in the strong 
optical pump of cavity mode a\. Considering that the fluctuations in the variables 
are all of the same order, and that Si always appears as a± + Si in the equations 
of motion for b and S2, we can ignore the dynamics of the pump fluctuations in the 
case where |oti| » (j&i), (j&z) and (b). This is the undepleted pump approximation. 
Adiabatically removing the pump from the dynamics of the system yields a pump- 
enhanced optomechanical coupling G = h\a.\\ between optical cavity mode a 2 and the 
mechanical resonance b. Dropping the subscript from the cavity mode a% and moving 
to a rotating reference frame results in the new effective Hamiltonian [To] . 

H eS = -Ab r b + Gatb + G*aP. (37) 

The system diagram and symbol corresponding to this simplified model of the PPT are 
shown in Figure |3j 
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Figure 3. Simplified PPT system diagram and symbol. The coupling rate G is 
proportional to h and \Jn\, where n\ is the number of photons in the pump mode. 



3.2. Scattering Matrix Formulation of the Phonon-Photon Translator 

To understand the properties of the PPT as a waveguide adapter, we begin with a study 
of its scattering matrix. Starting from the effective Hamiltonian given in equation (37), 



% The 7 in this case is in principle j e + 7^ as opposed to 7 e + 7i, and the equations must therefore be 
solved self-consistently. This is discussed in section 3.2.1 
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the Heisenberg-Langevin equations of motion for the Hamiltonian (37) are written under 
a Markov approximation in the frequency domain, 

— iua(u)) = — na(u) — iGb(u) — ^2K e a- m (u) — v / 2^a inii (u;), (38) 
— iu)b(uS) = — (7 — iA)b(u) — iG*a(uj) — -y/27 e 6 in (a;) — y / 2^& inji (u;).(39) 

The intrinsic noise terms a- mj i(u)) and &i nj i(u;) are the initial-state boson annihilation 
operators for the baths, while the extrinsic terms 5i n (o;) and bi n (u) are annihilation 
operators for the optical and mechanical guided modes for each respective waveguide. 



Since the effective Hamiltonian (37) has been used to derive equations (38 39), thus 



neglecting the counter-rotating terms present in the full system dynamics, the effects of 
phonon spontaneous emission noise is included separately. Following the discussion in 



Section 3.1, the effective Liouvillian (32) corresponds to replacing 7^ with 7^ in (39) and 



using a Langevin force &i n ,i(o;) satisfying the relations, 
(b- mA (ufbUu;)y=(n' + l)5(uj-uj% 



(40) 
(41) 



where 7^ and n' are given in equations (33 34). The intrinsic optical noise correlations 
are only due to vacuum fluctuations and given by (0,^(10') d\ n ; (o;)) = S(tu — to'). 

To ensure efficient translation, competing requirements of matching and strong 
coupling between waveguide and resonator must be satisfied. This is similar to the 
problem of designing integrated optical filters using resonators and waveguides [521 [53] . 
From the above equations and the input-output boundary condition [5"4"1 |4"9"] , we arrive 
at the matrix equation, 



flout (w) 
&out(w) 



s 



a- m (u) 
b- m (uj) 



flm,i(w) 



with scattering and noise matrices 



S 



su(u) s 12 (u) 
s 21 (u) s 22 (u) 



and Af 



n n (uj) u 12 (uj) 
n 21 (u) n 22 (uj) 



The elements of the scattering matrix S are 

2k s (j - i(A + u)) 



s u (uj) = 1 

Sl 2 (^) 

s 22 (u) 



\G\ 2 + ( 7 -z(o; + A))(k-iu)' 

2iG*^%K~ e 

\G\ 2 + (7 - i(u + A))(k - iuY 



'12 



2j e (K — icj) 



(42) 



(43) 



(44) 

(45) 
(46) 
(47) 



\G\ 2 + (7-1(0; + A))(k-zw)" 

Similar expressions are also found for the noise scattering matrix elements nij(u), with 
their extrema reported below. 

In order to obtain efficient conversion the cavities must be over-coupled to their 
respective waveguides, ensuring that the phonon (photon) has a higher chance of leaking 




Figure 4. Phonon-photon scattering matrix amplitudes for (7 e , 7$, K e , Ki, G) = 
2tt x (10, 1, 2000, 200, 155.6) MHz. (a) Plot of the frequency dependence of the optical 
reflection sn. The broad over-coupled optical line is visible, along with the PPT 
feature near the center in the unshaded region. This unshaded region is shown in more 
detail in plots (b), (c), (d), and (e), showing the frequency dependence of the scattering 

matrix elements Sn, S12, S21 an d S22, respectively. In each plot, the curves (— ), ( ) 

and (•— ) represent detunings of A = 0, 2ir x 200 and 2tt x 400 MHz respectively. 



into the waveguide continuum modes than escaping into other loss channels. In this 
regime, k % K e and 7 % j e . In the weak coupling regime, G < k, the response of 
the system exhibits a maximum for s 12 and S21 at u = and a minimum at the same 
point for Sn and 822- In fact, with realistic system parameters, only the weak-coupling 
regime leads to efficient translation. In strong-coupling (G > K and k » 7) the photon 
is converted to a phonon at rate G, and then back to a photon before it has a chance 
to leave through the much slower phononic loss channel at rate 7, causing there to be 
significant reflections and reduced conversion efficiency. To find the optimal value of G 
we consider the extrema given by 
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_ 2G^/'j e K e 



S 12 



P22 



G 2 + 7K 

G 2 + K7i - K'Je 



(49) 
(50) 

G 2 + K 7i + k% 

In the over-coupled approximation and in the case where Ki = ji = 0, it is easy to see 
that the full translation condition |si2| ma x = 1 is achievable by setting G equal to 

G = y/jK. (51) 

This result has a simple interpretation as a matching requirement. The photonic loss 
channel viewed from the phononic mode has a loss rate of — . Matching this to the 
purely mechanical loss rate of the same phononic mode, 7, one arrives at G° = <J^[k. 
The same argument can be used for the photonic mode, giving the same result. Under 
this matched condition, the linewidth of the translation peak in | S12 1 2 is simply 

4|G°| 2 

Ttransfer = = 2^. (52) 

K 

With intrinsic losses taken into account, either \su\ or | S22 1 (but not both) can 
be made exactly by setting G 2 = 7(k 6 — k^) or G 2 = (j e — respectively. The 
optimal state transfer condition, however, still occurs for G° = a/tk. The extremal 
values (u = 0) of the scattering matrix are in this case are, 

I I optimal ™i /E.o\ 

I I optimal _ le^e 

i^Uo -y— . ( 54 ) 

\s22UT 1 = (55) 

7e + 7i 

with corresponding noise matrix elements of, 



I I optimal _ \/ K i K e /c R \ 



I I optimal li^e ; r7 \ 

|ni 2 L P = o =A/— > (57) 

7e«t 
7K 



I I optimal jie^i / cc 

l"2i|«Ao =a/T7-' (5 J 



I .optimal _ V7»7e ,_„^ 
^22 w =0 ~ , ' l dy J 

For a set of parameters typical of an optomechanical crystal system, the magnitudes 
of the scattering matrix elements versus frequency are plotted in Figure |4j In these 
plots we have assumed resonant optical pumping of the a\ cavity mode and considered 
several different detuning values A. The normalized optical reflection spectrum (|sn| 2 ) 
is shown in Figure |2Fa), in which the broad optical cavity resonance can be seen along 
with a deeper, narrowband resonance that tunes with A. This narrowband resonance 
is highlighted in Figure |4](b), showing that the optical reflection is nearly completely 
eliminated on resonance. Photons on resonance, instead of being reflected, are being 
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converted into outgoing phonons as can be seen in resonance peak of |s2i| 2 shown in 
Figure [4](d) . A similar reflection dip and transmission peak is visible for the phononic 
reflection (S22) and phonon to photon translation (si 2 ) curves. It is also notable from 
Figure [4] that for small detunings A of the system (A < k, Cl) the resonant scattering 
matrix elements are only weakly affected and the translation process maintains its 
efficiency 



3.2.1. Modifications to Matching Condition due to Counter-Rotating Terms As a 
consequence of the counter-rotating terms treated in Section 3.1 , 7 is weakly dependent 
on G. In particular, by making the substitution 7^ — > 7-, where 7- is given by 

I G° P 

equation (33), the equation for the optimal G becomes = 7 e + 7^, where r ) i is 
itself dependent on G. Algebraic manipulations give us the desired value of G, 



G 



012 



simplifying in the sideband-resolved limit to G° 



(60) 



/7/c 



(l — (^) 2 ^ • The values for the 



scattering matrix elements given in equations (53 59|) are also suitably modified by the 
substitutions 7, —> 7^ and 7 —> 7 e + 7^ 



3.2.2. Phononic Waveguide Losses The issue of waveguide loss is one that is normally 
ignored in quantum optical systems where low-loss fiber or free-space links are readily 
available. On the phononic side of the systems studied here, the length of the waveguide 
and its intrinsic losses may be large enough for the waveguide attenuation factor to 
become important. Because of the negative effect of attenuation in applications involving 
optomechanical delay lines, and quantum state transfer, it is useful to model this loss 
and see how the scattering matrix elements are altered. 

A model of the system shown in Figure |5^b), where the lossy waveguide is replaced 
by a single beam splitter. This is accomplished by first modeling the lossy waveguide as 
a large number N of cascaded beam splitters, each reflecting 2aAz away from the main 
beam, with L = AzN. This serial array of beam splitters can be combined into one, 
with reflectivity 77 = e~ 2aL as iV —> 00. The relation for d out can be found by starting 
from the approximate scattering matrix relation for b out , 

Oout(t) = S 2 ia- m (t) + S22kn(t) + ni 2 kn,i{t) + 7l 2 2&in,i(*), (61) 

and using the beam splitter relations, 

d ou t(t) = V 7 / & out + V 1 _ V & wg,i,l (62) 
K(t) = VV d in + y/1 - 7] 6„g,i,2- (63) 

Assuming that the phonon cavity intrinsic loss bath is at the same temperature 
and uncorrelated with the phonon waveguide intrinsic loss baths, then one finds for 
translation through the lossy phonon waveguide, 

d ou t(t) = s' 21 d- m (t) + s' 22 d in (t) + n' 12 d- m;i (t) + n' 22 b inii (t) , (64) 
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where 



°21 

j 

S 22 



n 



12 



71 



Vv 821, (65) 

T) S 2 2, (66) 

Vv ni2, (67) 

V(l - v)v\ s 22\ 2 + ??|^22| 2 + 1 - V- (68) 

The value of S21 is simply reduced by a factor ^/r/ due to the lossy waveguide. For 
propagation lengths short relative to the attenuation length of the lossy waveguide, 
aL « 1, and the reduction in translation *Jrj % 1 — aL is small. The added noise due to 
the waveguide attenuation is contained in n' 22 , and is also seen to be small for aL « 1. 
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Figure 5. (a) System diagram of a PPT connected to a lossy waveguide of length L 
with attenuation factor a. (b) An equivalent system diagram in which the waveguide 
is replaced by a beam splitter with a reflectivity 77. 



3.3. Effects of Thermal and Quantum Noise 

In general, the type of noise which is relevant to the PPT depends on the conditions 
in which it is used. For example, when used as a bridge between RF-microwave 
photonics and optics for classical applications at room temperature, the thermal noise 
affecting the RF signal will be at a level which makes the quantum noise induced by 
spontaneous pump scattering irrelevant. On the other hand, when the system is used 
at cryogenic temperatures for connecting a superconducting circuit QED system to an 
optical system as described below, the quantum noise of the translation process itself 
becomes dominant. In what follows we analyze both sources of noise. 



3.3.1. Thermal Noise on the Optical Side The noise power on each of the output 
waveguide channels can be found by using the scattering matrix formulation described 



in Section 3.2 Here we evaluate the effects of thermal noise classically. When the 



system is used an optical drop filter, in which an optical beam is sent in and the optical 
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reflection is measured, we find: 

(t) = s n (t) * a in (t) + s 12 (t) * b in (t) 

+ n u (t) * a in ,i(t) + m 2 (t) * 6 in ,i(t) (69) 

where * represents convolution. Assuming that there are no cross-correlations between 
the various input noise terms, we find the classical spectral density of the noise to be: 

S'out(w) = \si 2 (u)\ 2 n + \n 12 (u)\ 2 (n + n spon ). (70) 

For a system with mechanical frequency less than 10 GHz, at room temperature 
(T = 300 K) the corresponding spontaneous emission noise is much less than the thermal 
noise, n spon « n ~ 10 3 , and n spon can be ignored. 



Ignoring the quantum noise term for the moment, and evaluating equation (70) 



at oj = after substituting in equations (56 59), the total thermal noise power on the 
reflected optical signal is found to be, 

P % = hu 2 n^nB, (71) 

K 

where B = 27 is the bandwidth of the PPT. In the classical regime, n = kT^^/^Ml), 
and this equation reduces to P^ ont = 2^^kT\, a ^{ui2/0){n e / n). This result has a simple 
interpretation. The noise power 7r7/cTb at h is the standard thermal noise input on the 
phononic side of the PPT. The ratio u 2 /Q is a translation factor, arising from the 
fact that quantum-limited conversion of phonons to photons causes an increase in 
energy by the factor of the ratio of their frequencies. Finally, the factor kJk is the 
extraction efficiency of photons from the optical mode to the waveguide. Alternatively, 
one may define an equivalent optical temperature by setting nn e kT oea = P^ out , yielding 
T ,eS = (2Qo/Qm)Tba&hi where the Q's represent the loaded quality factors of the optical 
and mechanical resonators. This last expression must be interpreted carefully, only in 
terms of a power equivalence, as spectrally the noise on the optical side of the PPT has a 
bandwidth of 27, while thermal noise radiating from the cavity would have a bandwidth 
of K. 

3.3.2. Phonon Spontaneous Emission Noise To calculate the effective increase in noise 
brought about by the spontaneous pump scattering and phonon emission process, we 



start from equation ( 70 ) and use the spontaneous emission contribution n spon to find, 




PL, = (^2)1»+ ; + i l(^)(2^-). (72! 

The second term in the brackets is due to the spontaneous emission of phonons by 
the optical pump beam, and is added to the thermal noise exiting the optical side of 
the PPT. Put in terms of an effective contribution to the thermal photon occupation 
number, the spontaneous pump scattering effectively adds, 

n °- spon = ~i (20) ' (73) 
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thermal photons to the cavity. This equivalence is only in terms of total noise power 
emitted, as spectrally the noise is emitted over the 27 bandwidth of the PPT resonance, 
not the entire optical cavity resonance as discussed above. 

4. Proposed On-Chip Implementation 

Up to this point in the analysis of the PPT, the discussion has been kept as 
general as possible. Such a system is, however, interesting only insofar as it is 
realizable, and we attempt here to establish the practicality of a PPT. Building upon 
recent experimental [17] and theoretical work [38l [39] , we provide the design of an 
optomechanical crystal formed in a silicon microchip that can realize a PPT system 
with high phonon-photon translation efficiency 

As previously mentioned, optomechanical crystals [T7J are engineered structures 
in which phonons and photons may be independently routed and their interactions 
controlled. In order to create a suitable OMC structure for the implementation of a 
PPT device one looks to a crystal lattice providing simultaneous phononic and photonic 
bandgaps for the guiding and co-localization of phonons and photons (55]. We have 
recently proposed [3H] such an OMC system, formed from a silicon-on-insulator wafer 
and consisting of a patterned thin membrane of silicon. The proposed "snowflake" 
crystal lattice supports a phononic bandgap in the 5-10 GHz mechanical frequency band 
and a photonic bandgap in the 1500 nm optical wavelength band. This quasi-2D crystal 
structure was also shown to support low-loss photonic and phononic waveguides, optical 
resonances with radiation-limited Q > 10 6 co-localized with mechanical resonances of 
frequency Q % 2tt x 10 GHz, and a single quanta optomechanical interaction rate of 
g % 2vr x 300 kHz. 

In this section we design an example OMC implementation of a PPT in the silicon 
snowflake crystal. We limit ourselves here to a two-dimensional (2D) crystal involving 
only 2D Maxwell's equations for transverse-electric (TE) polarized optical waves and 
in-plane elastic deformations of an infinitely thick slab. This simplifies the analysis 
and avoids some of the technical challenges related to achieving high optical Qs in 
a quasi-2D thin film structure, challenges which have already been studied and met 
elsewhere [3H EE1 EH] • 

4-1- Single and Double Cavity Systems 

The snowflake crystal consists of a hexagonal lattice of snowflake-shaped holes patterned 
into silicon as shown in Figure [6] The snowflake lattice used here is characterized by a 
lattice constant a = 400 nm, snowflake radius r = 168 nm, and width w = 60 nm. It 
possesses a full phononic bandgap from 8.6 to 12.6 GHz and a photonic pseudo-bandgap 
for TE optical waves from 180 to 230 THz. Defects in the crystal (features breaking the 
discrete translational symmetry of the underlying lattice) can support resonances with 
frequencies within the optical and mechanical bandgaps, leading to highly localized, 
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Figure 6. (a) The electric field magnitude |E| of the L2 cavity photonic resonance 
at v D = 199 THz. (b-e) The line defect optical waveguide structure, described in the 
text, which couples to the L2 cavity photonic resonance. The numerically simulated 
(c) acoustic bandstructure and (d) optical bandstructure of the line defect optical 
waveguide with W = 0.135a. (e) Plot of the out-of-plane component of the magnetic 
field, H z , of the guided optical mode at X-point of the bandstructure. (f) The 
displacement field |Q| of the L2 cavity phononic resonance at v m = 11.2 GHz. (g-j) The 
line defect acoustic waveguide structure, described in the text, which couples to the 
L2 cavity phononic resonance. The numerically simulated (h) acoustic bandstructure 
and (i) photonic bandstructure of the line defect acoustic waveguide with r c = 0.82r. 
(j) The magnitude of the mechanical displacement field, |Q|, of the guided acoustic 
mode at the X-point of the bandstructure. Calculations of the acoustic waveguide 
bandstructure are done using an FEM model [57], while for the optical waveguide 
simulations, a plane- wave-expansion method was utilized |58) . 
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strongly interacting resonances [39]. As an example, by filling an adjacent pair of 
snowflake-shaped holes in the crystal, a so-called "L2" cavity is formed which supports 
a localized photonic resonance at a frequency v a = 199 THz and a phononic resonance 
at a frequency v m = 11.2 GHz. The defect cavity structure and finite-element-method 
(FEM) [57] simulated field profiles of these co- localized modes are shown in Figures [6(a) 
and [6(f). 

The L2 cavity forms the basis of a more complex double-optical-mode cavity 
structure with the desired symmetry properties for efficient phonon-to-photon 
translation. By placing two separate L2 cavities close to one another, at sufficiently 
small separations, even and odd optical and mechanical super-modes form with split 
mechanical and optical resonant frequencies. We choose to displace the cavities from 
each other in the y direction, as shown in full system diagram of Figure [7] These super- 
modes of the coupled cavities are characterized with respect to their vector symmetry 
about the mirror symmetry transformation a y (x,y) = (x,—y). We denote these super- 
modes E+ and Q+, where ' + ' denotes symmetric and ' — ' denotes anti-symmetric 
symmetry with respect to a y . The supermodes can be written approximately as, 

E a + E b Q a + Q b 

E+ = -j= and Q+ = — -j= — , (74) 

where the subscripts a, b label the individual cavity fields. The cavity separation (14 
rows) shown in Figure [7] is chosen such that the optical supermode frequency splitting 
is very nearly identical to the mechanical mode frequency of v m = 11.2 GHz, as 
ascertained by FEM simulations. We focus here only on the mechanical mode of odd 
vector symmetry, Q_, since as we will show below, this is the mode which cross-couples 
the two optical super-modes to each other. 

4-2. Optomechanical Coupling Rates 

Optomechanical coupling (or acousto-optic scattering) arises from the coupling of optical 
cavity modes under deformations in the geometric structure. In the canonical form 
of radiation pressure, a mechanical deformation in the cavity induces a shift in the 
resonance frequency of a given cavity mode. The coefficient describing the cavity mode 
dispersion with mechanical displacement also quantifies the strength of the radiation 
pressure force that photons in the cavity mode exert back on the mechanical structure. 
More generally, mechanical deformations may couple one optical cavity mode to another. 

The self-coupling and inter-modal couplings caused by a mechanical deformation 
are modeled by the position dependent interaction rates gk(x) and h(x), respectively, 
in the Hamiltonian of eqn. pj). Both types of deformation-dependent optomechanical 
couplings may be calculated to first order using a variant of the Feynman-Hellmann 
perturbation theory, the Johnson perturbation theory [59J, which takes into account 
moving boundaries in electromagnetic cavities and has been used successfully in the 
past to model optomechanical crystal cavities [TTIEE]. The Hamiltonian, given to first- 
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order by, 



H = h^Uial&i + htttfb + - J]^,i( &t + 6 )«I%> ( 75 ) 



is then a generalization of that shown previously in equation (|2]), with 



-j-^dl (Q • n) (AeEf • E» - A(6^)D^ ■ Bj 



9ij = — \ / j (76) 

v mt y5^p|QIM^ e l E ilM^ e l E il 2 

where E, D, Q and t are the optical mode electric field, optical mode displacement 
field, mechanical mode displacement field, and the thickness of the slab or thin-film. 
For convenience, in what follows we denote the overlap integral in equation (76) as 
(Ej|Q|Ej). Since the optical supermodes are nearly degenerate, we replace Uij with 
either u + or w_ with little error, and following the notation used previously denote 
cross-modal coupling as h = g + _ = g~ t +. 

For the optical and mechanical modes of a single L2 cavity (see Fig. [6^a) and 
Fig. |6](f)) the optomechanical self-coupling is numerically calculated using ( [76] ) to 
be g = (E|Q|E)/27r = 489 kHz. Then using the symmetry selection rules in the 
overlap integrals, the only coupling terms involving the Q_ mode are the cross-coupling 
terms <E+|Q_|E T ). These terms are calculated to be <E+|Q_|E_) = <E_|Q_|E+) = 
(E|Q|E)/v2 to good approximation. For the supermodes of our double L2 cavity system 
this yields a cross-coupling rate of h/2ir = 346 kHz. 



4-3. Implementation of Waveguides 

A line-defect on an optomechanical crystal acts as a waveguide for light and sound |60, 
IBTj 1612] . In principle then, the same waveguide may be used to shuttle both the photons 
and phonons around on an OMC microchip. However, due to the different properties of 
optical and acoustic excitations, in particular, their typically disparate quality factors 
(roughly 10 6 for 200 THz photons and 10 4 for GHz phonons in silicon), the cavity loading 
requirements may be different for photons and phonons. As such, it is more convenient 
to implement two physically separate sets of waveguides, one for optics and the other 
for mechanics. 

Our chosen line defect for optical waveguiding, shown in Fig. [6^b), consists of a 
row of removed holes, with the rows above and below shifted toward one another by 
W such that the distance between the centers of the snowflakes across the line defect is 
\/3a — 2W. Simulations of this line defect waveguide show that there are no acoustic 
waveguide bands resonant with the localized mechanical L2 cavity modes of interest 
(see Fig. [6[c)), and therefore that this waveguide will not load the mechanical part of 
the L2 cavity system. Optically, this line defect has a single optical band crossing the 
frequencies of the localized optical L2 cavity modes of interest (see Fig. (6^d)), providing 
the single-mode optical waveguiding required for the PPT. 

The line defect used for the acoustic waveguide consists of a row of holes which have 
been reduced in size, as shown in Fig. [6Fg)). By shrinking the size of the central row 
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of snowflake holes by 18%, a single- mode acoustic waveguide is formed which spectrally 
overlaps the localized mechanical resonances of the L2 cavity (see Fig. [6|h)), and allows 
for mechanical coupling to the PPT. This same line defect waveguide is reflective for 
optical photons in the band of localized optical resonances of the L2 cavity (see Fig. 
(6fi)), and thus is isolated from the optical part of the PPT. Below we discuss how both 
the optical and mechanical waveguides may be used to load the PPT resonant cavity. 



4-4- Cavity-Waveguide Coupling 

By bringing the optical waveguide near the L2 cavity, the optical cavity resonance is 
evanescently coupled to the guided modes of the line-defect, as shown in Fig. [7^b). 
Control over this coupling rate is achieved at a coarse level by changing the distance 
(number of unit cells) between the cavity and waveguide. For the structure considered 
here, a gap of 8 rows is sufficient to achieve a coupling rate K e in the desired range. 
A fine tuning of the coupling rate is accomplished by adjusting the waveguide width 
parameter, with a value of W = 0.135a resulting in a loaded optical cavity Q- 
factor of Qwg,o ~ 3 x 10 5 (the corresponding external waveguide coupling rate is 
K e /2n = 300 MHz). Considering that intrinsic optical Q-factors as high as 3 x 10 6 
have been achieved in microfabricated thin-film silicon photonic crystal cavities similar 
to the sort studied here [63], the calculated optical waveguide loading should put such 
a cavity structure well into the over-coupled regime (n e / 'n^ % 10). A short section in 
which W is tapered is used to close off the waveguide on one side. 

The same design procedure for the acoustic waveguide results in an evanescently 
coupled waveguide at a distance of only one row from the L2 cavity. Since the acoustic 
line-defect waveguide does not support Bloch modes at the optical cavity frequency, 
no additional loss is calculated for the optical cavity resonance. In this geometry 
the mechanical cavity loading is calculated to be QwG,m = 1-3 x 10 3 , corresponding 
to an extrinsic coupling rate 7 e /27r = 4.4 MHz. Taking Qi % 10 4 as an achievable 
intrinsic mechanical Q-factor (we need a citation to the Berkeley group here for Si GHz 
resonators), such a loading also puts the mechanical system in the over-coupling regime 
(7e/7, * 10). 

Simulations of the above cavity-waveguide couplings are performed using FEM [57] 
with absorbing boundaries at the ends of the waveguide. The resulting time-average 
electromagnetic Poynting vector |(S )J 2 = |E x H*|/2 of the optical field leaking from 
the L2 optical cavity resonance is plotted in Fig. [7](b), while the mechanical Poynting 
vector |(S m )J 2 = | — v-T| 2 (v is the velocity field, and T the stress tensor) of the acoustic 
waves radiating from the mechanical mode of the L2 cavity is shown in Fig. [7tc). It is 
readily apparent from these two plots that the coupling of the two different waveguides 
to the L2 cavity act as desired; the acoustic radiation is coupled only to the phononic 
waveguide, and the optical radiation is coupled only to the photonic waveguide. 

In order to individually address and out-couple from the even and odd symmetry 
cavity resonances of double-L2-cavity structure used to form the PPT, a pair of 
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Figure 7. (a) The full optomechanical crystal PPT system, consisting of a pair 
of coupled L2 defect cavities with acoustic and optical waveguide couplers. The 
waveguide coupling to the cavities, for both optics and acoustics, consists of a pair of 
horizontal line defect waveguides, one to each of the L2 cavities. The optical waveguides 
(highlighted in blue) are the two outer waveguides and the acoustic waveguides 
(highlighted in beige) are the two inner waveguides. The relevant symmetric and 
anti-symmetric optical (blue) and acoustic (beige) supermodes of the cavity- waveguide 
system are shown as envelope functions with the appropriate symmetry (dashed for 
symmetric and solid for anti-symmetric). The FEM-simulated (b) photonic (|(S ) t | 2 ) 
and (c) phononic (|<S m ) t | 2 ) Poynting vectors of the lower cavity-waveguide structure, 
illustrating the selective optical loading of the lower waveguide and the selective 
mechanical loading of the upper waveguide on the cavity. 
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waveguides is used for each of the optical and mechanical couplings. As shown in 
the overall PPT design of Figure [7|a), each of the L2 cavities has an acoustic and 
an optical waveguide coupled to them. Excitation of a pair of waveguides either in 
or out of phase would thus allow for coupling to the symmetric or anti-symmetric 
supermodes, respectively, of the double-L2-cavity. Similarly, spatial filtering (via an 
integrated directional coupler or waveguide filter for instance) of the output of a pair 
of waveguides would allow for the selective read-out of either the symmetric or anti- 
symmetric cavity modes. One could in principle utilize spectral filtering to perform 
the selective mode coupling; however, with the narrowband nature of the optical and 
mechanical supermode splittings, spatially independent channels of excitation and read- 
out may be a preferred option. 

In summary, the OMC PPT as designed couples the symmetric and anti-symmetric 
optical modes of a double-L2-cavity system via a co-localized anti-symmetric mechanical 
resonance at frequency v m = 11.2 GHz. In the notation of Section [3j the lower frequency 
symmetric optical mode is the pump mode (cavity mode ai), the anti-symmetric mode 
is the signal mode or cavity mode 02, and the anti-symmetric mechanical resonance 
is phonon mode b. Both optical modes are designed to have a resonant frequency 
in the near-IR around a frequency of 200 THz, with a frequency splitting engineered 
to be equal to the mechanical frequency, w_ — u + % Q = 2n x 11.2 GHz. The 
numerically calculated waveguide and optomechanical coupling rates for this system 
are (ft e ,7 e ,/i) = 2tt x (300,4.4,0.35) MHz, with the required number of intracavity 
pump photons for optimum operation {G % ^Jl^rf e ) of such a PPT estimated to be only 
|c!i,o| 2 = 1-1 x 10 4 (assuming minimal intrinsic losses and n e « k, 7 e ~ 7). 

5. Applications 

At the simplest level, the extremely narrow optical response of the PPT, as shown 
in Figure |4j provides the opportunity for design and fabrication of filters with MHz- 
scale linewidths in the optical domain. By comparison, a purely passive optical design 
would require optical cavities with quality factors of Q % 10 8 . More generally, such 
a scheme demonstrates a promising aspect of optomechanics in the realm of optical 
information processing. In this section three examples applications of the PPT are 
studied in detail. The first two, optical delay lines and wavelength converters are 
further examples of optical information processing that is of considerable interest in 
both classical and quantum information processing. The last application, using the 
PPT to provide optical "flying qubit" capability to superconducting microwave quantum 
systems, is an example of how optomechanics can have a fundamental role in hybrid 
quantum system engineering. 
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Figure 8. Schematic layout (left) and signal path (right) of a PPT used in an optical 
delay line/filter configuration in which a phonon mirror at the end of the acoustic 
output reflects outgoing phonons back into the PPT. The large optical delay is afforded 
by the acoustic path length of the signal in which acoustic waves propagate some 10 5 
times slower than photons. Optical filtering is provided by the narrow resonance 
bandwidth of the mechanical component of the PPT. 



5.1. Delay Lines 

Efficient reversible conversion between traveling photons and phonons can be used to 
realize an optical delay line, as shown schematically in Figure [8j In this geometry, 
the acoustic waveguide used to extract phonons from the PPT is terminated abruptly, 
forming an effective acoustic wave mirror which reflects phonons back toward the PPT 
after some propagation distance and delay Resonant photons sent into the optical port 
of the PPT, will then re-emerge, reflected and delayed by twice the length of the acoustic 
waveguide. The delay line functionality comes from the inherent slowness of acoustic 
waves in comparison to electromagnetic waves (roughly a factor of 10 5 for waves in 
silicon). For a similar reason, electro-acoustic piezoelectric structures are used to create 
the chip-scale RF-microwave filters found in many compact wireless communication 
devices |64|. 

The usefulness and bounds on the main characteristics of a PPT-based delay line, 
i.e., the total delay possible and the delay-bandwidth product, may be simply estimated 
without referring to a particular implementation of the system. The maximum possible 
delay is given by the lifetime of an excitation on the mechanical side of the system 
(cavity and waveguide), and is given by I/tj, which is limited by material properties of 



the mechanical system. The bandwidth of PPT conversion is given by eqn. (52), and is 
the total loss rate seen by the mechanical resonance, 27. As such the delay-bandwidth 
product can at most be 

Aut ~ — , (77) 

li 

which is approximately twice the acoustic waveguide to mechanical resonance over- 
coupling ratio 7 e /7i in the PPT. Manipulation of the acoustic waves within the delay 
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waveguide itself, before conversion back to optics, may also be envisioned, enabling 
existing phononic information processing capabilities [33J to be applied to optical signals. 



5.2. Wavelength Conversion 




Figure 9. (a) Schematic of a back-to-back PPT structure for optical wavelength 
convertsion/delay. For pure wavelength conversion, the phononic waveguide is 
vestigial, and can be removed by simply coupling the optical elements in both PPTs 
to the same mechanical resonant clement. The simplified effective system diagram for 
such a device is shown in (b), where optical resonances a and c are coupled to the same 
mechanical resonance b. The pump cavities for each system are omitted. 



Figure [9j(a) shows the schematic for a PPT based device that can act simultaneously 
as a narrowband filter, a delay line and a wavelength converter. It consists of connecting 
serially by a common acoustic waveguide, two PPT devices operating at different optical 
but matched mechanical frequencies. Interestingly, if the only goal is to perform photon- 
to-photon wavelength conversion, one can omit the connecting acoustic waveguide. By 
placing the two PPTs "on top of each other", such that the optical cavities on both 
PPTs are coupled to the same mechanical resonance, as shown in Figure [9]^b), optical 
wavelength conversion can be accomplished. Such a PPT geometry could be realized 
by either having photonic cavities with multiple modes, or using two photonic cavities 
coupled to the same mechanical mode. Such photon-to-photon conversion could even be 
taken to an extreme, allowing, for instance, optical-to-microwave wavelength conversion 
if one of the photonic cavities is a microwave cavity. 

For the simplified wavelength conversion system of Fig. P(b), the PPT matching 



condition of equation (51) and noise analysis of Section 3.3 carry over with only minor 
adjustments. For the simplified wavelength conversion system the thermal noise is now 
split between the two optical channels (a and c of Fig. [9jb)), while the spontaneous 
emission noise in the system is approximately doubled (for similar optical cavities) 
due to the two uncorrelated spontaneous emission processes occurring from the optical 
pumping of each individual cavity. In a single element PPT, the optimal G matches the 
pure mechanical damping of the mechanical resonance (7) to the induced optomechanical 
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loading of the mechanical resonance (G 2 / k) by the optical cavity. The matching 
condition for the simplified wavelength converter now must balance a mechanical 
resonance coupled on one side to an optical cavity with induced mechanical loading 
rate G 2 Jn a and on the other side to a second optical cavity with induced loading rate 
G 2 /n c . As such, assuming that 73 « G\/k^ we arrive at the photon-photon converter 
matching condition, 

K a K c 

The optomechanical system as described would act as a quantum-limited optomechan- 
ical wavelength converter. Finally, we note that this particular implementation of the 
wavelength converter could also function in a wider array of optomechanical platforms 
since there is no longer a need for phononic waveguides. 

5.3. Quantum State Transfer and Networking between Circuit QED and Optics 
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Figure 10. System diagram for quantum state transfer from optical to 
superconducting qubits. The two 3-port components are circulators with each input 
connected to the output directly clockwise from it. 



Two of the key requirements for a viable platform for quantum computation 
are the ability to do storage and communication of quantum information. For the 
case of superconducting phase qubits, promising theoretical proposals to provide 
such functionality have involved electrical [65] and mechanical resonators (H |5]. 
Experimentally, an electromagnetic resonator quantum bus was demonstrated by 



Proposal for an Optomechanical Traveling Wave Phonon-Photon Translator 



26 



Sillanpaa et al. [66] in 2007, while more recently O'Connell et al. [ID] demonstrated 
the strong coupling of a mechanical resonator to a superconducting qubit. Circuit 
QED (cQED) to date remains limited by the lack of a true long range state transfer 
mechanism, one which is readily available for the case of quantum-optical qubits, in the 
form of optical fibers and free-space links. Using the PPT system, one potentially could 
implement a version of the quantum state transfer protocol of Cirac et al. [67] allowing 
for the high fidelity transfer of states between optical and superconducting qubits. Such 
a system would satisfy one of the original goals of hybrid quantum system [68] by 
interfacing a quantum optical and solid-state qubit. 

By connecting the phononic waveguide of a PPT to a piezoelectric resonator 
strongly coupled to a superconducting qubit [U El [10], and connecting the PPT at its 
optical end to an optical cavity QED system, as shown in Figure 10 , the phonon-photon 
translator could be used as an intermediary in a state transfer protocol among two 
energy-disparate quantum systems. The optical system (A) is composed of a Fabry- 
Perot cavity containing a three-level atom system in a A configuration. As shown 
in [67], in the correct limit, this system may be modeled with effective Jaynes-Cummings 
(JC) Hamiltonian with a Rabi frequency controlled by another beam. For the 

superconducting system (£>), a mechanical resonance is coupled to a phase qubit, with 
a bias current used to change the frequency of the resultant two-level system, which 
effectively changes the coupling rate between the qubit and mechanical resonator. This 
leads again to a system with an externally controllable Rabi frequency of gsif)- The 
Hamiltonian of each subsystem is then given by 

Hj = hg^e-^ajd] + hg^e+^afa (j =A,B), (79) 

where dj are the annihilation operators of the photonic or phononic resonances external 
to the PPT, and frj are the level lowering operators for the respective qubits to which 
they are coupled. 

Each cavity mode, with annihilation operators Cj, is coupled to its respective 
waveguide with a loss rate Fj. To characterize the PPT, the intrinsic losses in these 
cavities are ignored. Additionally, the phononic waveguide is assumed to be loss- 
less, though losses may be taken into account through the minor readjustment of 
the scattering parameters studied in Section |3.2.2| Using the input-output boundary 
conditions \5Q the frequency domain expression for the noise input into systems A 
and B are found to be, 

CA, in(w) = ain,c(w), (80) 

CB, m(w) = S 2 l(uj)^/f~AC A {uj) + S a i(w)ai njC (w) + S 22 M&in,c M 

+ n 12 (u)ain,i(uj) + n 22 (u)b inti (uj), (81) 

where a injC and 6 irijC represent the noise being coupled into the system from the third 
input of each circulator in Figure 10 To convert this equation to the time-domain, the 
convolution between various operators and scattering matrix elements must be taken. If 
the photon pulse is of sufficiently large temporal width, i.e., with bandwidth less than the 
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bandwidth of the PPT, the frequency dependence of each scattering matrix element can 
be removed, replacing it with its extremal value assuming that the system is operating 
at resonance (A = 0). This requires that the coupling rates gj(t) should change slowly 
relative to the response of the PPT. Under this condition, the input-output relations in 
the time-domain are then 

CA, in(t) = a irL7C (t), (82) 

CB, in(t) = S21\/TAC A (t) + S 2 l6in,c(0 

+ s 2 2&in,cW + n i2a- m ,i(t) + n 2 2&in,i(£)- (83) 

In modeling the noise of the system, it is assumed that the optical noise inputs are 
in the vacuum state, and that the phononic noise is thermal with thermal phonon 
occupation numbers n and n' for &m iC (i) and b\ n ^(t) , respectively, where the PPT phonon 
spontaneous emission noise is combined with the intrinsic thermal bath coupling of the 
mechanical mode in n' as described in section 13. II 

Using standard operational methods of Quantum Stochastic Differential Equa- 
tions [HI EH], the master equation describing the evolution of systems A and B is 
found to be, 

P= \\- H a + h b, p] + ^^i,aP + ~y£i,bP 



+ ~y (\ s 22\ 2 n + |n 22 | V) C t ,bP 

+ \/TaTb\s 2 i\ ([4, c A p] + [pc\ } c B ]) (84) 
with the Liouvillians Cij and Ctj are given by, 

Cijp = 2cjpc] - c]cjp - pcJjCj, (85) 
C T jP = Ci,aP + 2cjpcj - Cjc]p - pcjc]. (86) 



The final term in the master equation ( 84 ) is the cascading term [SU EHJ EO] which gives 
rise to the unidirectional coupling between the systems. 

For the PPT, parameters typical to an OMC structure such as the one 
with scattering matrices plotted in Fig. [4] are used, (j e , ji, n e , Q) = 2n x 
(10, 1, 2000, 200, 8000) MHz. Equation can be used to find the optimal matching 
optomechanical coupling rate, which for the assumed PPT parameters is G° = 
155.4 MHz. At this optomechanical coupling rate, the resonant noise and scattering 
matrix parameters of the PPT are (s 2i , s 22 , n 21 , n 22 ) = (0.917,0.074,0.290,0.262), with 
a spontaneous emission noise equivalent occupation number of n spon = 0.200. Assuming 
that the PPT is cooled to the same cryogenic temperature of the superconducting qubit 
system (T < 100 mK), the thermal bath component of the effective thermal occupancy 
of the PPT mechanical resonance can be neglected, and n' ~ 0.251. 

The exact functional form of the gj(t) are found through numerical optimization. 
This was done by taking the pulse-shape to be a step smoothed by a sinusoidal function 
with a rise (fall) time of t r (tf). Optimization on the state transfer fidelity for an 



ideal PPT (s 21 = 1) and for circulators running in the direction shown in Fig. 10 
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(qubit transfer from optical to superconducting system) leads to rise and fall times of 
tA = 23.5 fis and = 16 [is, which are within the modeled PPT's (27)/27r = 22 MHz 
bandwidth. 




Figure 11. The main plot shows the time variation of the fidelity F\ found by 
numerically solving the quantum master equation for a version of the quantum 
state transfer protocol introduced by Cirac, et al. [67] for the state transfer from 
the optical to superconducting qubit. PPT parameters of (■y e , K e , Ki, Q) = 2tt x 
(10,1,2000,200,8000) MHz are used here to find the the optical (F^ A ; — ) and 
superconducting (Ff :B ; — •) state-transfer fidelity curves for an initial state of |1)a|0)b. 

A B 

We see that a maximum fidelity F 1 ' of approximately 0.8 for the final state transfer 
is possible. The continuous lines (— ) and (— ) were found by quantum simulation of an 
ideal PPT, i.e. 7, = Kj = 0. In all cases analyzed the external system parameters were 
(T Al T B ,g A la *,g™*) = 2ir x (50,5.0,5.0,1.0) MHz. The inset plot shows the shape of 
the control pulse used for the optical ( ) and superconducting ( — ) qubits. 



Putting this all together, in Fig. [IT] we plot estimates of the fidelity of the quantum 
state transfer between system A and B via the connecting PPT. The definition of 
fidelity used to calculate the state transfer efficiency is F^' B = r R: a(\j b)(j b\p) ■, where 
j = 0, 1, + represent respectively the ground |0), excited |1), and |+) = 2~ 1//2 (|0) + |1)) 
states of the atomic and superconducting two-level systems. Under these conditions, and 
considering an ideal PPT (s 21 = 1), states are transfered with fidelities F 1 ,B = 0.9983, 
F+' B = 0.9995 and Fq' b = 1.00. Taking into account the actual scattering and noise 
matrix values given above for the PPT, and accounting for the spontaneous emission 
noise of the PPT we find that the fidelities are reduced to F x ' = 0.803, F^ a = 0.936 
and Fq' d = 0.983. The inverse system, with circulators turning the opposite direction 
to transfer qubits from the superconducting to optical system was also studied, for which 
the same input pulses only time reversed, and yield fidelities F B ' A = 0.772, F B ' A = 0.904 
and F BA = 0.983. 
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6. Summary 

We have introduced the concept and design of a traveling phonon-photon translator. 
We have shown that with a realistic set of parameters and the use of existing silicon 
optomechanical crystal technology, efficient and reversible conversion between phonons 
and photons should be possible. By characterizing the noise processes experienced by 
such a device, both classically and quantum mechanically, we have shown the utility 
of traveling phonon-photon translation to important problems in both classical optical 
communication and quantum information processing. 
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